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Abstract. iKhomenko et al.l estimate the mean magnetic field strength of the quiet Sun to be 20 G. The figure is smaller than 
several existing estimates, and it comes from the comparison between observed Zeeman polarization signals and synthetic 
signals from numerical simulations of magneto-convection. The numerical simulations require an artificially large magnetic 
diffusivity, which smears out magnetic structures smaller than the grid scale. Assuming a turbulent cascade for the unresolved 
artificially smeared magnetic fields, we find that their unsigned magnetic flux is at least as important as that explicitly shown 
in the simulation. The unresolved fi elds do not produce Zeeman polarization but contribute to the unsigned flux. Since they are 
not considered bv IKhomenko et all their mean magnetic field strength has to be regarded as a lower limit. This kind of bias is 
not specific of a particular numerical simulation or a spectral line. It is to be expected when observed quiet Sun Zeeman signals 
are compared with synthetic signals from simulations. 
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1. Rationale 

The observational properties of the qui et Sun magnetic 
fields are not well established (see, e.g., ISanchez Almeidal 
12004 and references therein). We are in the way of char¬ 
acterizing them, a process whose end is difficult to predict 
partly due to contradictions among existing measurements. 
In this sense, IKhomenko et al.l ll2005ll estimate a mean field 
strength or unsigned flux density {B) of only 20 G, a value 
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20061). This Note argues that the unsigned flux estimated by 
Khomenko et alJ should be regarded as a lower limit rather 
than an unbiased estima te. If the conclusi on is correct, it may 
reconcile the result of IKhomenko et akl with other existing 
estimates claiming {B) to be of the order of 100 G. Moreover, 
the argumentation and the caveat it entails are of general 
nature. They must be considered whenever the quiet Sun 
unsigned flux is derived by comparison of observed Zeeman 
signals with numerical simulations. Such approach to study 
the quiet Sun magnetism is bound to play a major role in the 

future. _ 

The paper bv IKhomenko et all ll200-5l) compares the polar¬ 
ization observed in a quiet Sun region with the synthetic po- 


Send offprint requests to: J. Sanchez Almeida, e-mail: jos@iac. es 


larization produced by realistic numerical simulations of solar 
magneto-convection. The numerical simulations have a mean 
magnetic field strength {B) that decreases exponentially with 
time f, 

<B(f)> - <B(0)) exp(-f/TB), (1) 

with an e-folding time scale td - 50 min (Khomenko & Vogler 
2005, private communication). As usual, the angle brackets 
denote volume average. The snapshot of the simulation with 
{Bit)) ^ 20 G is the one that best reproduces the histograms 
of observed polarization signals. This agreement is used to 
indicate that the observed inter-network regi on has (B) close 
to 20 G. There is a natural way to reconcile IKhomenko et alJ 
il200.5t) result with the much larger fluxes from other measure¬ 
ments. The numerical simulations do not provide information 
on structures smaller than a few grid points. They miss small 
scale magnetic structures that would show up if the simula¬ 
tions had finner more realistic resolution. The simulations only 
include a part of the magnetic fields, and so, it is to be expected 
that the estimate by IKhomenko et alJ ll2005l) only provides a 
lower limit to the true solar {B). This consideration can be a 
serious caveat or an academic remark depending on the im¬ 
portance of the neglected magnetic fields. There is no definite 
way of evaluating the importance of the bias except for improv¬ 
ing the spatio-temporal resolution of the simulations to reach 
realistic values. This approach, however, demands a computer 
power exceeding by orders of magnitude the current resources. 
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Our work provides a preliminary estimate using the tools avail¬ 
able at present. 


2. Energy and unsigned flux density in the 
unresoived magnetic fieids 


The simulations assume a magnetic diffusivity rj, - 1.1 x 
10 ^' cm^ s“', which is some three orders of magnitude larger 
than the Ohmic diffusivity to be expected at the base o f the pho¬ 
tosphere (77 ^ 10® cm^ s”'; iKovitva & Cramlll983h . This un¬ 
realistically high diffusivity is needed for technical reasons to 
stabilize the numerical MHD code, and to prevent the artificial 
built-up of energy at the grid s cale ( IVogleil2003HVogler et al.l 
boost l^rdlund & Steinll99^ . How can the numerical simula¬ 
tions be realistic if they employ such artificially large magnetic 
diffusivity? They can if there is a physical mechanism that dif¬ 
fuses away the magnetic energy at the rate imposed by the arti¬ 
ficial diffusivity. Somehow the existence of such mechanism is 
implicit when the numerical simulations are used to represent 
and characterize the real Sun. Lacking of obvious alternatives, 
we will describe such mechanism as turbulent diffusion, sim¬ 
ilar to that required for the astrop hysical dynam os to operate 
within reasonable time scales (e.g., IParkeJ 19791 Chapter 17). 
This association allows us to carry out simple estimates. The 
turbulent diffusion is produced by a complex unresolyed mag¬ 
netic field with spatial scales / so small that the diffusion time 
scale for the true small Ohmic diffusiyity is similar to the dif¬ 
fusion time scale for the structures in the simulation jParkeJ 
fl97i lBiskamp<2003l) . Assuming that the high magnetic diffu¬ 
siyity used by the numerical code is of turbulent nature, one can 
eyaluate the magnetic energy in the turbulent magnetic field im¬ 
plicit in the simulation. This complex magnetic field does not 
produce polarization signals, so that the obse ryables remain as 
those synthesized by iKhomenko et akl il200.'ili . Neyertheless, it 
contains energy and unsigned flux that must be included when 
estimating the true unsigned flux density corresponding to the 
obseryed polarization signals. The section aims at showing how 
such turbulent magnetic field can indeed contain a significant 
amount of energy and unsigned flux. 

In order to estimate the magnetic energy in the turbu¬ 
lent cascade from the minimum resolyable scale L, to the 
diffusiye scale /, we need a model for the MHD turbulence. 
The MHD turbulence is a topic of actiye research and there 
is n o unique and final way to approach the problem (see, 
e.g., lRoldvre^ll20^■5^ [Brandenburg et ah boo.^l and references 
therein). Howeyer, one can estimate the magnetic energy in the 
turbulent cascade using a few approximate prescriptions which 
are now in use. They all lead to the conclusion that this energy 
may be significant. 

We start off by assumi ng a Kolmogoro y spectrum for 
the turbulent cascade te.g.. iBiskamnlEfif)^ . Then the total 
magnetic energy per unit mass corresponding to a (spatial) 
wayenumber k is. 




( 2 ) 


The symbol e stands for the energy dissipated per unit time and 
unit mass by the turbulent cascade (see, e.g., lBiskamr)ll2003l 
§ 5.3.2). The Kolmogoroy constant Ck is a numerical factor 
of the order of 1.7. Then the magnetic energy per unit mass in 
unresolyed fields £, is, 

E,^ (4) 

JiniL 2(2nYi^ 

where we haye taken into account that I L. Now, the mean 
magnetic field of the numerical simulations haye an exponen¬ 
tial decrease with an e-folding time scale To (equation IQ). 
Assuming that the standard deyiation among the field strengths 
scales with the mean field*, 

(5) 


then. 


{B^} - cr| + {Bf oc {B)\ 


(6) 


implying a magnetic energy {B^)f(^7T) decreasing exponen¬ 
tially with a time scale half td (compare equations Q and 0 ). 
The energy dissipated per unit time by the turbulent cascade 
turns out to be 


E 


p dA Stt ^ 


-1 -I 


(7) 


with p representing the mass density required to transform en¬ 
ergy per unit mass {Et and EY) into energy per unit yolume 
{B^)l{^n). The dissipated energy has a compact expression in 
terms of the Alfyen speed va, 


..2 _-i 


(8) 


with. 


47Tp 


(9) 


If the symbol Bt stands for the turbulent magnetic field strength, 
(B?) _ <b2) v2 


Szrp <B2) 2 ■ 

Using equations 0 , and (O, 


( 10 ) 


^ 3Ck / L \ 2/3 
<B2) (2;r)2/3lTB va/ 


( 11 ) 


The yariables defining our problem are, p ^ 3 x 10^’ g cm ® 
(the mass density at the base of the photosphere), = 

28 G (assuming ctb — {B} ^ 20 G), L ^ 100 km (AnnendixiAl. 
and To - 50 min (equation 0). They render. 


<B 2 ) 


0.56. 


( 12 ) 


We are interested in the mean magnetic field strength (or un¬ 
signed flux density), since this parameter is used to characterize 


with 

27t/L < k < 2nll. 


' A condition typical of many probability density functions suitable 
for describing the variable B, e.g., a Maxwellian distribution, and ex- 
(3) ponential distribution, a uniform distribution, etc. 
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the observations. Assuming the relationship 0 with the same 
ratio o-bUB) for both {B) and (B,), 


{B) 



(13) 


The analysis above has been repeated using the so-called 
Iroshn ikov-Krishnan spectrum of MHD turbulence dRiskamnl 

[ 2 ^, 

Ek^C,K(vAe)^'^k-^'^- (14) 


(MHD turbulence simulations seem to yield s pectra in betwee n 
Kolmogorov and Iroshnikov-Krishnan; see iBoldvre^ bOOSt l 
Following the procedure explained above, one finds, 


(Bj) ^ 2 Qk/ l y/2 


(15) 


With CiK - Ck, and for the numerical values used in equa¬ 
tion <12> . 


<B 2 > 


1.30, 


(16) 


yielding. 


{B) 


^ 1.14. 


(17) 


The estimates above assume the inertial range of the turbulence 
to begin at a scale larger than L. For this reason we can describe 
the magnetic energy for k > InjL using laws characteristic of 
the inertial range (equations 0 and ini). However, this as¬ 
sumption is not decisive. Consider a turbulent cascade whose 
inertial range begins at an unresolved scale A < L. Let us rep¬ 
resent the unknown spectrum in the range between L and A as 
a power law of index a. 


Ek 


{B^),L.y-a 

Snp 


r“. 


(18) 


with 2nlL < k < InjA. The factor multiplying kr°‘ guaran¬ 
tees the continuity of Ek at k - InjE, it is equal to the mean 
energy per unit mass and spatial frequency existing in the re¬ 
solved scales, so that for k — InjL, 


(B^) 1 

Sttp InjE 


(19) 


The exponent a is also imposed by continuity arguments. When 
the inertial range begins, i.e. at k = 27r/A, then Ek is given by 
equation (E). This constraint leads to 

(a - 1) ln(L/A) = In (^) - In C, - ^ In (^). (20) 

The dehnition of £, in equation 0 implies, 

pin IK pl-^l^ pin IA 

E, - j Ekdk + j Ekdk > j Ekdk ^ 

*JlnjL *Jln/A UlnjL 


(B^) 1 

Snp a — I’ 


( 21 ) 


where we have assumed a > 1 and L » A. (Intermediate val¬ 
ues of A provide magnetic energies in between equation ED 
and the limit considered here.) The expression (ED, together 
with equation (Gni», render 



1 . 


( 22 ) 


Once more the unsigned flux in unresolved fields turns out to 
be significant. The ratio in equation (I22t has been evaluated 
using equation (ED with the same parameters leading to equa¬ 
tion (ED, which provide an exponent a varying from 2.1 to 1.8 
when A goes from 10 km to 0.1 km. 

Equations (I13> . (I17t and \22l suggest that the implicit tur¬ 
bulent held accounting for the used magnetic diffusivity has a 
magnetic flux similar to that explicitly shown in the simula¬ 
tions. 

Some numerical simulations of turbulent magneto¬ 
convection show a range of wavenumbers where the magnetic 
ener gy exceeds the k i netic energy (super e quipartition range; 
see iBiskamnl l2003t iMaron et^ l2004t iBrandenburg et alJ 
l2005h . The magnetic energy densities used aboye are much 
smaller than the kinetic energy density of the solar granula¬ 
tion. If this super equipartition range would haye to be included 
in our estimate, one needs to increase e with respect to the 
yalues employed aboye. The turbulent magnetic energy would 
increase accordingly, leading to {B,) larger than the yalues in 
equations ED and ED- Consequently, the estimates aboye are 
probably conseryatiye, a conclusion reinforcing the importance 
of the unresolyed magnetic fields. 


3. Conclusions and discussion 

The numerical simulations of solar magneto-conyection require 
artificially large magnetic diffusion, which smears out all the 
spectrum of magnetic structures smaller than the grid scale. 
The use of such high magnetic diffusiyity can be understood 
as the effect of a complex unresolyed turbulent magnetic held 
with spatial scales so small that the diffusion time scale for 
the true small Ohmic diffusiyity is similar to the diffusion time 
scale for the structures in the simulation. Such an implicit mag¬ 
netic held does not contribu t e to th e polarimetric signals syn¬ 
thesized by lKhomenko et ai]ll2005l) . Assuming a turbulent cas¬ 
cade for the unresolyed artificially smeared magnetic fields, we 
And that their unsigned magnetic flux is at least as important 
as that explicitly shown in the simulation. Should this mag¬ 
netic flu x is considered, th e Zeeman polarization signals mea¬ 
sured by lKhomenko et aD are consistent with an unsigned flux 
of 2 X 20 G (i.e., (B,) ^ (B); equations 11.31. 1171 and 1221 V 
In othe r words the unsigned flux assigned by lKhomenko et aH 
(I 2 OO 5 I) must be regarded as a conservative lower limit. This 
conclusion is not specific of the simulations analyzed here. A 
bias is to be expected whenever the quiet Sun unsigned mag¬ 
netic flux is inferred as the unsigned flux of numerical simula¬ 
tions reproducing observed Zeeman polarization signals. 

The calculations that we describe represent only a first ap¬ 
proximation to estimating the bias. They are based on the the¬ 
ory of MHD turbulence, which remains to be completed. The 
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conclusions have to be backed up or rejected by numerical sim¬ 
ulations with realistic Ohmic diffusivities. 
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Appendix A: L - 100 km 

The sp atial discretization of the MHD equations solved by IVdgleJ 
120031) is based on a four-order finite difference scheme. Then the 
components of the magnetic field are approximated by a four-order 
polynomial, so that the first and second partial derivatives appearing 
in the MHD equations are represented by a third-order polynomial 
and a second-order polynomial, respectively. We will estimate L as 
the smallest wavelength of a sinusoidal magnetic field for which the 
error involved in these polynomial approximations is insignificant. 

Let us denote by f{x) one of the components of the magnetic field 
vector, with k the wavenumber along the spatial coordinate x. 


f{x) = sin(fcx). 


(A.1) 


The numerical code represents it in the surroundings of xo as the poly¬ 
nomial f„(x). 


./«W = £,/'*"*(xo) 


(x - Xo)" 


(A.2) 


The fifth order term in the Taylor expansion of /(x) provides the error 
of the approximation, A/(x), 


= /(x) -/a(x) . (A.3) 

As usual, the symbol /*"' denotes the «-th derivative of /(x). The ap¬ 
proximation is carried out at each pixel, so that 


|x - Xol < Ax, (A.4) 

with Ax the pixel size. Among the various polynomials used to repre¬ 
sent the MHD equations, the computation of the second derivatives, 

/®(x) =-fc^/■(x), (A.5) 

corresponds to the lowest order and, therefore, it is the worst approx¬ 
imation. If the error associated with the computation of the second 
derivatives is tolerable, then the polynomial approximation is tolera¬ 
ble. According to equation the error Ap^\x) is 

A/*^’(x) = f^\x) - jf\x) k^ cos(kxo) (A.6) 

where we have taken into account that /*^'(xo) = cos(kxo). The 
average (unsigned) error within the range of the approximation (equa¬ 
tion lA.TI I is, 

_ 1 pA'o+A-r 

|A/ra(x)i = — |A/®(x)|rfx 

2Ax 

I cos(kxo)|—-. (A.7) 


The factor | cos(kxo)| varies with Xq but its average is the same as the 
average of |/(x)|. Then the mean relative error when computing the 
second derivatives, 6, is given by, 

_ ^3 A „3 

|A/ra(x)l/l/‘'>(x)| = 5 - (A.8) 

This expression provides a relationship between the relative error of 
the approximation, the wavelength of the sinusoidal A, and the pixel 
size, 

27T 

A = — ^4.7 Ax (6/0. 1)^‘A. (A.9) 

k 

We define L as the smallest wavelength for which the error is insignif¬ 
icant, namely, A assuring 6 < 0.1. According to equation and 

keeping in mind that Ax = 21 km, 

L 4.7 Ax 100 km. (A. 10) 
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